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ABSTRACT 

Physics of the gravitational effect of the galactic bar and spiral structure is pre- 
sented. Physical equations differ from the conventionally used equations. 

Application to the motion of the Sun is treated. The speed of the Sun is taken to 
be consistent with the Oort constants. 

Galactic radial migration of the Sun is less than ±0.4 kpc for the four-armed 
spiral structure. The Sun remains about 75 % of its existence within galactocentric 
distances (7.8 — 8.2) kpc and the results are practically independent on the spiral 
structure strength. Thus, the radial distance changes only within 5% from the value 
of 8 kpc. 

Galactic radial migration of the Sun is less than ± (0.3 — 1.2) kpc, for the two- 
armed spiral structure. The Sun remains (29 — 95)% of its existence within galactocen- 
tric distances (7.8 — 8.2) kpc and the results strongly depend on the spiral structure 
strength and the angular speed of the spiral arms. The radial distance changes within 
(3.8 - 15.0)% from the value of 8 kpc. 

If observational arguments prefer relevant radial migration of the Sun, then the 
Milky Way is characterized by the two-arm spiral structure. 

Key words: Galaxy - galactic bar - galactic spirals - motion of the Sun - equation 
of motion. 



1 INTRODUCTION 

Orbital evolution of the Sun in our galaxy, Milky Way, is im- 
portant for understanding of galactic evolution, its dynam- 
ics, kinematics and chemical composition. It is also relevant 
for understanding of evolution of the Solar System, e.g. the 
effect of the galactic tides on the evolution of the Opik - 
Oort cloud of comets. Several papers on the gravitational 
effect of the galactic bar on motion of the Sun in the Galaxy 
appeared recently, e.g. Minchev & Famaey (2010), Famaey 
& Minchev (2010), Minchev et al. (2010), Minchev et al. 
(2011). One of the important results of the papers shows 
that the Sun can radially migrate in the disk of the Galaxy. 

Minchev & Famaey (2010), Famaey & Minchev (2010) 
present the change of the solar galactocentric distance. The 
distance may change within several kiloparsecs. 

If we are interested in real galactic radial migration of 
the Sun, we have to treat physical equation of motion. The 
first aim of this paper is to find the equation of motion, to 
put the relevant equations on a firm physical basis. The sec- 
ond aim of this paper is to use the physical equations in per- 
forming calculations similar to Minchev & Famaey (2010). 
We are interested in radial motion (radial migration) of the 
Sun in the Galaxy. 

Section 2 presents multipole expansion of the galactic 



potential and the concentration is paid to the quadrupole 
term, important for the galactic bar. The relevant quantities 
found in the monopole and quadrupole terms are discussed 
in Section 3, where the values of the quantities are calculated 
for several models of the galactic bar. Section 4 offers the 
equation of motion for an object under the gravitational 
influence of the galactic disk, halo and the monopole and 
quadrupole terms of the galactic bar (the monopole term 
corresponds to the conventional idea of the galactic bulge). 
Section 5 presents numerical results for the radial migration 
of the Sun. Finally, discussion and conclusion can be found 
in Sections 6 and 7. 



2 MULTIPOLE EXPANSION 

Having a mass density distribution g(r), the gravitational 
potential at a point r is 

m = -G J |^d 3 r', (1) 

where G is the gravitational constant. Making an expansion 
in terms of the ratio = r'/r, we can write 
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and, 



$(°)(r) 
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g(f') d 3 



G M 



(2) 



(3) 



where M is the total mass concentrated in a small region so 
that the condition \f '\/\r\ <C 1 is fulfilled. Other terms of 
the expansion are (f) = 0, and, the quadrupole term is 
(Einstein's summation convention is used) 

G n 

2 r 6 



$ (2) (r) 



[3(r')i (r')j - \r'\ 2 S t j] g(f') d 3 r' , 

n. = (r)i/r (4) 

and 5ij is the Kronecker delta. Since & 2 \r) is given by a 
contraction of two tensors of the second order, Qij and rii rij , 
it is an invariant independent on the coordinate basis. The 
numerical calculations can be done in a special coordinate 
basis, the primed coordinates, e.g., defined by the relations 
describing rotation 

x = + x cos a + y sin a , 
y = — x sin a + y cos a , 

z = z . (5) 

In the case of a galactic bar we will use 

a = Q. b t + a , (6) 

where \ilb\ is the magnitude of the angular velocity of the 
bar's rotation. In reality, fij, is negative and the negative 
sign denotes the observational fact that rotation is negative, 
i.e., clockwise. 

Eqs. (J)-© yield 

G 

<*> (2) (x,y,z) = - — -3 x X# , 

2 r 5 

X<s> = Qi\n\ + Q 2 2nl - (Qn + Q22) n\ , (7) 

where primes above the Qn and Q22 terms are omitted, for 
the purpose of brevity. The coordinate system is chosen in 
the way that Qy = if i 7^ j and some symmetry exists. The 
mass density is an even function of coordinate arguments, 



p(x, y, z) (2x 2 — y 2 — z 2 ) dxdydz , 
Q22 = / p(x, y, z) (2y 2 — x 2 — z 2 ) dxdydz (8) 



and 



2 
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2 
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(2 2 ,2-2 l0 • 

la; cos a + y sm a + Ixy sm a cos a) fr , 

/ 2 . 2 , 2 2 o • "\/2 

Ix sm a + y cos a — Ixy sm a cos a) fr , 
z 2 /r 2 . (9) 



Inserting Eqs. (£9|) into Eq. ([7]) one obtains 

^ (x, y ,z) = _gqu(i-q M /qu) 



4 r° 

„ I + Q22/Q1I / 2 . 2 o 2\ 

+ (x 2 — j/ 2 ) cos (2q) + 2 2; y sin (2a) , 
a = Qb t + Qo , (10) 

if also Eq. © is added. As for the Milky Way, the value of ao 
is about —25 degrees (compare with Majaess 2010; the real 
value may differ in more than 10 degrees) and the value of 
the angular velocity is fib = — 55.5 km s _1 fcpc -1 (compare 
with Minchev & Famaey 2010), t is the time. 



2.1 Discussion on Eq. (|10p 

Our result represented by Eq. (|1U|) differs from the result 
conventionally used by other authors. The conventional re- 
sult is (see, e.g., Minchev & Famaey 2010, Dehnen 2000) 



*ioL = | Qt v 2 cos [2(0 - fl b t)] X X<I,2 C , 



1 



X$2c 
Xij>2c = 2 



m 



r^r b 



(5) ■ 



r < 7-f, 



(11) 



where r b = 3.44 fcpc, u c = 240 fcm s _1 , 0.1 < Qt < 0.4 
(Minchev & Famaey 2010). 

At first, the second formula of Eqs. (j 1 1 p does not hold. 
The inner part of the galactic bar potential must be calcu- 
lated from the Poisson equation or its equivalence: A$ = 
4tvGq, $ = — G J g(f ')/\r — r'\ d 3 r '. As for the case 
r ^ rf, in Eqs. (|11[) . it leads to 

3 . 



®co„v(x, V,z) = - Q T 



2 \VR 2 + z 2 J R 2 

Y&2c = {x 2 — y 2 ) cos(2tlbt) 

+ 2 x y sin(2fM) . (12) 



We remind that r = V R 2 + z 2 , where 7? is the 2-dimensional 
radial coordinate in the galactic equatorial plane and z is 
the vertical coordinate normal to the equatorial plane. One 
can immediately see that physical Eq. (j 1 Of) differs from the 
conventional potential. 



3 MODELS OF THE GALACTIC BAR 

As examples of mass distribution in the galactic bar we take 
the models presented by Stanek et al. (1997). For the mod- 
els Gl, G3, E2, PI, P2, P3 we introduce the following com- 
plete nonsymmetric ellipsoid coordinates (r is a dimension- 
less quantity, now - it differs from the quantity used in the 
previous section): 

x — r xq sin 9 cos cp , 
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y = r yo sin 6/ sm<p , 
z — r zq cos 8 , 

with the corresponding volume element 
dx dy dz = xo Vo zo r 2 sin d6 dip dr . 



(13) 



(14) 



The corresponding quantities Qn and Q22, and, also mass 
of the galactic bar M, can now be calculated in an analytical 
way. 

(i) The model Gl with mass density 



PGi(x,y, z) = p exp I 



72), 



(-Y+ - 

\xo/ \ yo 



z 
zo 

yields 

M = 2 V2 tt 3/2 p x y zo , 
Qn = 2 V2 7r 3/2 po 2:0 Vo z (2x1 - Vo - z o) 

= M (2x1 - Vo - Zo) j 
Q22 _ 2yl -x\- Zq 
Q11 2x 2 -Vo- z o ' 

(ii) The model G3 defined by the mass density 
pG 3 (x,y,z) = po r-" 1,8 exp(- r 3 ) , 
produces 



(15) 



(16) 



(17) 



M = 4-k x 0.7394 x p x yo z , 

4tt / 2 2 2\ 

Qn = — x 0.3219 x po xo yo zo (2x - y - z ) 
= 0.1451 x M {2x1 



2 2\ 

Vo- Zo) , 



Q22 _ 2yl -£§ - zl 
Qu 



2xq J/ 2 Zq 



(18) 



where the number 0.3219 is an approximate value of the 
integral 



2.2 . 
r exp I 



1 dr . 



(iii) The model E2 defined by the mass density 
PE2(x, y, z) = po exp (- r) 
gives 

8tt po xq yo z , 



(19) 



M 
3u 



327r po xo yo z (2xq - yo - zl) 



2/Q11 = 



= 4 M (2xo 
2y 2 



- xl - 



2 

- Vo 
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2\ 
Z ) 
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(iv) The model PI with mass density 
Po 



p P1 (x,y,z) 



(1 + rf 



(20) 



(21) 



yields 



M = 47T po x yo zo 



(1 + vY 



■ dr 



In 



r 3 

— po xo yo zo bar 3 , 
6 (1 + r bar ) 



4t /„ 2 2 2^ 

— po xo yo zo [2x - y - z J 



(1 + rY 



■ dr 



4tt 



Po 



xo yo zo (2xq - yl - zl) x X P1 



M (2x1 - Vo - zo) 



3r bar + 22 + 



30 

Tbar 



+ 



12 



12 



fl + — ) ln(l + r 6ar ) 

V Thar J 



Thar 



3(l + r bar ) 3 
rb/xo (as an example) . 



- 4 In (1 + Tfcar) 

(22) 



We have taken the symbol r bar as the upper limit of the 
integral in Eq. 1)220 (if one would like to use infinity as the 
upper limit of the integral, then the integral would diverge, 
Q11 — > 00. If we take into account that the real length of the 
bar r b is proportional to xo, then we can take r bar = r b /xo- 
(v) The model P2 with mass density 



1 \ Po 
p P2 (x,y,z) = 

r (1 + r) 

yields 



M = 47r po xo yo zo 



= 27T po x yo zo 



4tt 



(23) 



(1 + rY 



■ dr 



r 2 

' bar 

(1 + r bar 



— po x yo zo (2x'l - Vo - zo) 



4tt 



{1+rY 



■ dr 



j Po x yo zo (2x1 - Vo - zo) x X P2 

6 



- M (2xt - yl- zl 



2r bar + 9 + 



r ba 



6(l + —) \n(l + r bar ) 



X f 



2r Lr + jHLr + £> T bar 
2 (1 + rbarf 



3 In (1 + rb ar ) 



r bar = r b /xo (as an example) . 
(vi) The model P3 with mass density 
Po 



p P3 (x,y,z) 



(1 + r 2 )^ 



(24) 
(25) 
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yields 

M = 47T po xo 2/o Zo 



dr 



Qi 



(1 + r 2 ) 2 
o 

2ix po x y z ^arctan r bar - 

47T / 2 2 2\ 

— po xoyozo [2x - J/,, - z Q ) 



^ po x yo z (2x1 - 2/o - z o) 



7*bar 
1 + r lar 



(1 + r 2 ) 2 



1 Tbar 



2 1 + r£, p 



■ arctan rt, ar 



3 M (2x1 



2 
J/0 



2 

Zo) x 

rtar + r bar / [2 (l + r 2 or )] - (3/2) arctan r bar 
arctan r bar - r bar / (1 + r 2 ar ) 

rt ar = r b /xo (as an example) . (26) 

As for the model El with mass density 

PEi(x,y,z) = p exp(-r e ) , (27) 
where 

re = M + M + M, (28) 

xo yo zo 

the values of M, Qn a Q22 can be easily analytically calcu- 
lated (one may avoid Egs. I13l - I14[) : 

M = 8 po x yo z , 
Q11 = 16 p x yozo (2x1 - Vo - z o) 
= 2M (2xl-yl-zl) , 

o 2 2 2 
2j/o - Xo - Zo 



2/Qll = 



The model G2 with mass density 
p G 2(x,y,z) = p exp (~r 2 s /2) , 
where 



\a;o/ 



+ 



1/1 



+ 



(29) 
(30) 

(31) 



can be solved with the following substitutions: 

• 1/2 a 

X = Xo fs sin 17 COS </J , 

■ 1/2 n ■ 

J/ = 1/0 fj sin f sin (p , 

1/2 a 

z = zo r s cos 17 , 

for the upper hemisphere, and, 

■ 1/2 a 

x = xo r s sin ' v cos ip , 

■ 1/2 q ■ 

y = jo r, sin 17 sin ip , 



(32) 



Z = - zo r s cos 



1/2 







(33) 



for the lower hemisphere, and, for both hemispheres r s 6 
(0,oo), 1? G (0, 7r/2), <£> £ (0, 2tt), and, the volume element is 

dx dy dz = — a;o 2/o 20 r 2 cos _1//2 9 dO dip dr s (34) 

for both hemispheres. The values of the quantities M, Qn 
and Q22/Q11 are: 



M 
Qn 



Q22 
Q11 



2 7r 3/2 if f-^j po xo yo z , 
po x yo z x Xg2 



2y 2 ~x 2 - x 



M x X G 2 , 
2^[r(3/4)] 2 Zo 2 



2«o -2/o - V 2 /tt [r(3/4)] 2 z 2 



X G 2 = 2x Q - y 



(I) 



2 

2() 



(35) 



where K(x) is the complete elliptic integral of the first kind. 
Finally, the model E3 



pE-i(x, y, z) = po Ko(r s ) , 
yields 



(36) 



M 



po xo yo z 



18 n po x yo z x X E 3 

-1NT-1 



9V2 



K 



G) 



Q22 

Qn 



22/0 ~ xl 



M x X B 3 , 
2A[r(3/4)] 2 z 2 



2x1 -y 2 - VV" [ r (3/4)]^o 2 



o 2 2 

2z ~ 2/0 



(!) 



2 



(37) 



3.1 Numerical results 

Table 1 presents numerical results for the data given by 
Stanek et al. (1997 - Table 4). Since the integrations in 
the models PI, P2 and P3 lead to integrals which require 
finite limits in order to be convergent, we have considered 
the upper limit corresponding to the length of the bar r b 
and the values of xo for individual models (see Eqs. 1221 1241 
I26[) . Taking into account the value r b — 3.44 kpc, the re- 
quired quantities are presented in Table 2. If we admit a 
relative error of r b to be 20 %, then the relative errors of the 
quantities are also presented in Table 2. 



4 CONVENTIONAL AND PHYSICAL 
APPROACHES 

In order to compare the published results with our approach, 
we have to summarize important equations. 
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Table 1. Numerical values of the quantities discussed in the pre- 
vious subsection. Data for xq, yo and zq are taken from Stanek 
et al. (1997). The quantity M/po contains the central density po. 





Model 


X() 




7n 






[pc] 


[pc] 


[pc] 




Gl 


1478 


754 


611 




G3 


4537 


2303 


1725 




E2« 


1000 


500 


400 




E2< 2 ' 


1000 


505 


388 




PI 


1565 


783 


582 




P2 


2653 


1320 


975 




P3 


1869 


939 


713 




El 


1650 


868 


474 




G2 


1298 


694 


572 




E3 


1075 


563 


442 












Model 


Q22/Q11 


Q 


hi/M 


M/po 




[-] 




pc 2 ] 


[pc 3 ] 


Gl 


-0.4146 


3.4271 x 10 6 


1.0724 x 10 10 


G3 


-0.3938 


4.7727 x 10 6 




E2« 


-0.4151 


6.3600 x 10 6 


5.0266 x 10 9 


E2( 2 ) 


-0.4017 


6.3777 x 10 6 


4.9245 x 10 9 


PI 


-0.3957 








P2 


-0.3957 








P3 


-0.3999 








El 


-0.3224 


8.9338 x 10 6 


5.4309 x 10 9 


G2 


-0.4461 


5.7115 x 10 6 


1.0639 x 10 10 


E3 


-0.4294 


1.2084 x 10 7 


6.9228 x 10 9 



Table 2. Numerical values of quantities for the models PI, P2 and 
P3. The values of Qu/M and M/po hold for r b = 3.44 kpc. The 
relative errors q(Qh/M) and g(M/po) hold for the uncertainty 
of 20 % in r b . 



Model 


Q22/Q11 


Qu/M 


Q{Qu/M) 




[-] 


[pc 2 ] 


[%] 


PI 


-0.3957 


2.3149 x 10 6 


-31 / + 34 


P2 


-0.3957 


2.1628 x 10 6 


-32 / + 36 


P3 


-0.3999 


2.5152 x 10 6 


-29 / + 31 




Model 


Q22/Q11 


M/po 


s(M/po) 




[-] 


[PC 3 ] 


[%] 


PI 


-0.3957 


0.9700 x 10 9 


-20 / + 17 


P2 


-0.3957 


6.8383 x 10 9 


-19 / + 16 


P3 


-0.3999 


5.1389 x 10 9 


-22 / + 18 



4.1 Conventional approach 

The conventional approach to galactic motion of the Sun 
defines the following potentials of the Galaxy (Minchev et 
al. 2010). The background axisymmetric potential due to 
the disk and halo has the form 



corresponding to a fiat rotation curve. The non- 
axisymmetric potential due to the galactic bar is modeled 
as a pure quadrupole 



(2) 



- Q T vl cos [2(0 - n b t)] x X (2) 

* b cc 



(3 



3" 



, r < n 



(39) 



where r b = 3.44 kpc, v c = 240 km s" 1 , 0.1 < Qt < 0.4 
(see also Minchev & Famaey 2010) and the spiral potential 
is given by 



cos [a ln(r/j-rj) — m(cf) — Q s t)] 



(40) 



where e s is the spiral structure strength, 0.015 < |e s | < 0.072 
for a two-armed spiral structure and 0.007 < |e s | < 0.031 for 
a four-armed spiral structure, e s — sign{— a) x |e s |, a = —4 
and —8 for the two-armed and four-armed spiral structure, 
respectively, where the negative sign corresponds to trailing 
spirals, m is an integer corresponding to the number of arms, 
n s G {0.7,0.8,0.9, 1.0, 1.1} x v c /R , i.e., Q s [km s" 1 kpc' 1 ] 
e {21, 24, 27, 30, 33}, if v c = 240 km s' 1 and R = 8.0 kpc. 



4.2 Physical approach 

In this section we take the background axisymmetric po- 
tential due to the disk and halo in the form similar to Eq. 
(1381) 



4> 



disk-\-halo 



v c {1 - m(r/r G )} 



(41) 



where ra is the dark matter radius of the Galaxy. Eq. (|4ip 
is the physical potential corresponding to the flat rotation 
curve defined by the condition v(r) — v c . The result pre- 
sented by Minchev & Famaey (2010) or Dehnen (2000), see 
our Eq. ([38]). differs from Eq. (|4T|) . 

The potential due to the galactic bar is given by Eqs. 
© and CES): 



<$>bar(f) 

$ ((,) (r) 
<1> (2) (x,j/,z) 



$ (0) (r) + $ (2) (f) , 

J g(f') d 3 f' = 
G Qn (1 - Q22/Q11) 



G M 



4 r 5 

1 + Q22/Q11 



I 2 , 2 o 2 
(x + y — 2 z 



4>o 



v c ln(r/r ) 



(38) 



M 

hi[M@ pc 2 ] 
Q22/Q11 
a 
ao 
6 



I-Q22/Q11 
+ (x 2 — y 2 ) cos (2a) + 2xy sin (2a)] , 
1.4 x 1O 1O M , 
8 x 10 6 x M[M Q ] , 

- 2/5 , 

fl b t + ao , 

- 25° , 

- 55.5 km s _1 kpc -1 , (42) 
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where also the value of M from Dauphole & Colin (1995 - 
Table 5) is used and the values from Table 1 are taken into 
account . 

The spiral potential is given by the conventional model 
represented by Eq. (|40[) . in a slightly corrected form (Rb = 
3.44 kpc): 



s corr 



— e s Vc cos [a In ( R/Rb 



<j>sO + O s t)] 



v\ | Si cos (jn(fi a o + & In — 
+ 5*2 sin (rnfiso + a In -^-^j | 



Si 
S 2 

cos (40) 
sin (40) 

cos (20) 
sin (20) 



cos (m<f>) cos (mf2 s i) — sin (m0) sin (mtt s t) 
sin (m0) cos (mfi s i) + cos (m0) sin (mQ s t) 

4 



(I)' + (I)' 



= 4 



a; 


y 






i? 


R 


w - 


(1 



(l) 1 - 



2*1 



sO = «0 , 



(43) 



where the quantity s o denotes orientation of spiral arms 
in the coordinates x-y at the intial time t = and we put 
4> a o — Qo- The quantity e s is the spiral structure strength, 
0.015 < |e s | < 0.072 for a two-armed spiral structure and 
0.007 < |e s | < 0.031 for a four-armed spiral structure, e s 

■ x |e s | , a — + 4 and + 8 for the two-armed and four- 
armed spiral structure, respectively, where the positive sign 
corresponds to trailing spirals, m is an integer corresponding 
to the number of arms, Q a £ {0.7, 0.8, 0.9, 1.0, 1.1} x v y o/Ro, 
where v y o is the initial circular velocity of the Sun and Ro 
= 8.0 kpc. 

4.3 Equation of motion 

If an object moves in the Galaxy under the action of galac- 
tic disk, halo, galactic bar and spiral structure, then the 
equation of motion of the object reads 

v = - grad {^disk+haio + &b ar + $s corr) , (44) 
where Eqs. (liTj) . (|42)l and J43j can be used. 

4.4 Motion of the Sun — initial conditions 

Initial conditions for the Sun are 

x = Ro — 8 kpc , 

V = o , 
v x = , 

v y = v y o = — 255.2 km s _1 , 



Table 3. Numerical values of minimum and maximum galacto- 
centric distances of the Sun under the action of galactic bar, disk 
and halo. Motion for the time span of 5 X 10 9 years is considered. 





Vy 




Rmax 


[ km s _1 ] 


[ km s _1 ] 


[ kpc ] 


[ kpc ] 


240.000 


-255.196 


7.384 


9.193 


202.176 


-220.000 


7.774 


8.317 



for Eqs. gU O - the first case , 
v y = v y o = — 220.0 km s _1 , 

for Eqs. ((41] [42l [43]) — the second case 
v y = v y o = + 240.0 km s _1 , 

for Eqs. (EU [39] [40]) , 



(45) 



E(0 = 



(K^)] 2 ) s(0) 



- V^ 2 + S(w 2 ) 

(b(^)] 2 ) $(0) + (b(%)] 2 ) 4)( 



*(2) ' 



i? 



-Ro 

G Qn 



GM 
"RrT 



1 + 



_3_ 

32 



1 + ^1 x 



4i? 3 ' On 
GM 



Ro ' 



Ro 
(46) 



for the physical approach. Eqs. (|38[) -(145 [ ) uniquely determine 
motion of the Sun in the Galaxy. 



5 COMPUTATIONAL RESULTS 

We numerically solved equation of motion in accordance 
with Eq. (|44l) and the results presented in the previous sec- 
tion. We used both of the physical values of v y , as it is given 
by Eqs. (1451 - the first and the second cases). 



5.1 Effect of the galactic bar 

Table 3 shows the extremal values of the evolution of the 
galactocentric distance of the Sun under the action of galac- 
tic halo, disk and galactic bar. The results presented in Table 
3 hold for the special case when the effect of spiral struc- 
ture is ignored. While the nonexistence of the nonsymmetric 
parts in term $^ 2 ' would produce circular orbit for the ini- 
tial conditions given by Eqs. (|45l - the first and the second 
cases) , R — Ro — 8 kpc, the nonsymmetric parts in the term 
$ ( - 2 ^ cause slight oscillations in R(t). 

The relative change in R is about 15 % (or even larger) 
according to Famaey & Minchev (2010), Minchev & Famaey 
(2010). Our physical approach yields two different values. 
The first case, when v y o = —255.196 km s _1 , the initial 
orbital velocity of the Sun, yields Rmin/Ro = 0.923 and 
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Table 4. Numerical values of minimum and maximum galacto- 
ccntric distances of the Sun under the action of galactic bar, disk 
and halo. Motion for the time span of 5 X 10 9 years is considered, 
Vy0 = -220.0 km s _1 . 



Qu/M 






Rmax 


[PC 2 ] 


[ km s~ 1 ] 


[ kpc ] 


[ kpc ] 


3.4271 x 10 6 


202.742 


7.752 


8.300 


1.2048 x 10 7 


201.669 


7.796 


8.333 



Rmax/Ro = 1.149. While the value of Rmax corresponds to 
the value given by Famaey & Minchev (2010), the value of 
Rmin is greater than the value presented by other authors. 
Moreover, we have to stress that the effect of spiral arms con- 
sidered by Minchev & Famaey (2010), Famaey & Minchev 
(2010), Minchev et al. (2010), Minchev et al. (2011) is not 
considered in our calculations. The second case, when v y o 
= -220.0 km s _1 , yields R min /R = 0.972 and R max /Ro 
= 1.040. The galactocentric distance of the Sun changes in 
much smaller interval than it was presented in the previous 
papers. Again, the effect of spiral arms is not considered in 
our calculations. 

Let us consider motion of the Sun during 5 x 10 9 years. 
The first case, characterized by the value v c — 240 km s~\ 
yields that Sun moves within galactocentric distance R G 
(7.75 kpc, 8.25 kpc) during 44.3% of the total integration 
time, and, R G (7.5 kpc, 8.5 kpc) during 71.5% of the to- 
tal time. The second case, characterized by v y o = —220.0 
km s _1 , yields that Sun moves within galactocentric dis- 
tance R G (7.9 kpc, 8.1 kpc) during 50.5% of the total time, 
and, R G (7.8 kpc, 8.2 kpc) during 83.2% of the total time. 
The results do not change if we consider the total time span 
10 x 10 9 years instead of 5 x 10 9 years. 

In reality, there exists some dispersion in the values of 
Qn. Taking into account the results of Table 1, we can con- 
sider the following two extremal values: Qu/M = 3.4271 x 
10 6 pc 2 for the model Gl, and, Qu/M = 1.2048 x 10 7 pc 2 
for the model E3. The results for the radial migration of the 
Sun are summarized in Table 4. The Table 4 shows that the 
values are: Rmin ^ 7.75 kpc and Rmax ^ 8.33 kpc. 

5.2 Simultaneous effect of the galactic bar and 
the spiral structure 

If we take into account Eqs. (|43p . then galactic radial mi- 
gration of the Sun is less relevant than Minchev & Famaey 
(2010), Famaey & Minchev (2010), Minchev et al. (2010), 
Minchev et al. (2011) present. This holds both for the two- 
and four-armed spiral structures, although the effect of the 
two-armed spiral structure is much more relevant than the 
effect of the four-armed spiral structure, see Tables 5 and 7. 

The results summarized in Tables 6 and 7 show that 
radial migration of the Sun is practically unaffected by the 
effect of the four-armed spiral structure with v c — 202.176 
km s _1 . The results are practically independent on the spiral 



Table 5. Numerical values of minimum and maximum galacto- 
centric distances of the Sun under the action of galactic bar, disk 
+ halo and spiral arms. The two-armed spiral structure and v c 
= 202.176 km s -1 , v y0 = -220.0 km s" 1 and v x0 = 0.0 km s" 1 
are considered. The results hold for the time integration of 5 X 
10 9 years. The Sun remains (T / 100) X 5 X 10 9 years in the 
galactocentric distances Ft G (7.8, 8.2) kpc. 









Rmax 


T 


[-] 


[ Vyo/Ro ] 


[ kpc ] 


[ kpc ] 


[%] 




U. i 


/ .oyo 


O.ZDD 


qa n 




U . o 


/. /oy 




no q 

yo.o 


U.Ulo 








(1Q 1 

yo.i 




1.0 


7.758 


8.232 


94.5 




1.1 


7.744 


8.224 


92.8 




0.7 


7.434 


8.445 


52.7 




0.8 


7.632 


8.248 


76.5 


0.030 


0.9 


7.673 


8.197 


82.6 




1.0 


7.658 


8.141 


74.9 




1.1 


7.696 


8.164 


85.4 




0.7 


7.109 


8.752 


39.0 




0.8 


7.434 


8.362 


51.1 


0.045 


0.9 


7.519 


8.298 


60.4 




1.0 


7.509 


8.225 


59.3 




1.1 


7.582 


8.200 


67.5 




0.7 


6.772 


8.971 


29.3 




0.8 


7.201 


8.534 


37.0 


0.060 


0.9 


7.339 


8.459 


44.2 




1.0 


7.333 


8.346 


44.3 




1.1 


7.437 


8.268 


50.0 



structure strength. The found interval for radial migration 
of the Sun (7.8 — 8.4) kpc is much shorter than the result 
presented by Minchev & Famaey (2010), Famaey & Minchev 
(2010), Minchev et al. (2010), Minchev et al. (2011). We can 
explain the difference by incorrectness of the equations used 
by the above mentioned authors. Moreover, it seems to be 
nonphysical to use the case of resonance Q, s — v y o /Ro as the 
relevant for the motion of the Sun in the Galaxy. There is 
no argument why the spiral arms should rotate with the an- 
gular velocity equal to the angular velocity of the Sun. The 
well-known result states that a flat rotation curve exists for 
the Galaxy. The flat rotation curve yields an orbital speed 
practically independent on the galactocentric distance and 
the angular velocity is given by the ratio of the orbital speed 
and the galactocentric distance. Thus, the angular speed of 
a galactic object is a function of the galactocentric distance 
and the fixed value of Q s , Q s = v y o/Ro, is a pure mathe- 
matical toy and it contains no relevant physics. 

The last column in Tables 5-7 denotes the relative time 
T of the existence of the Sun in the galactocentric distances 
R G (7.8, 8.2) kpc. The two-armed spiral structure yields 
T < 50% for large values of spiral structure strength. The 
four-armed spiral structure yields T G (69%, 79%). 
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Table 6. Numerical values of minimum and maximum galacto- 
centric distances of the Sun under the action of galactic bar, disk 
+ halo and spiral arms. The four-armed spiral structure, v c = 
202.176 km s' 1 , v y o = —220.0 km s~ 1 and v x q = 0.0 km s~ 1 
are considered. The results hold for the time integration of 1 X 
10 9 years. The Sun remains ('/' / 100) X 1 X 10 9 years in the 
galactocentric distances Ft S (7.8, 8.2) kpc. 





U s 




Ftmax 


T 


[-] 


[ Vyo/R0 } 


[ kpc ] 


[ kpc ] 


[%] 




U. / 


i . /oy 


q on7 
o.zy / 


77 ft 
1 /.b 




U . o 


7 77/1 

/. / /4 


o.oll 


7Q ft 
/o.b 


U.UU / 


u.y 


*7 7ft n 

/. /oy 


o.oUo 


70 c 
(O.O 




1.0 


7.779 


8.300 


78.2 




1.1 


7.772 


8.310 


78.8 




0.7 


7.764 


8.299 


72.4 




0.8 


7.769 


8.328 


75.4 


0.014 


0.9 


7.766 


8.322 


74.7 




1.0 


7.783 


8.309 


74.8 




1.1 


7.770 


8.326 


75.8 




0.7 


7.758 


8.329 


69.9 




0.8 


7.765 


8.345 


73.9 


0.021 


0.9 


7.763 


8.334 


74.2 




1.0 


7.781 


8.323 


73.2 




1.1 


7.768 


8.339 


74.3 




0.7 


7.752 


8.360 


69.1 




0.8 


7.759 


8.381 


72.7 


0.028 


0.9 


7.761 


8.357 


72.4 




1.0 


7.778 


8.340 


72.9 




1.1 


7.763 


8.355 


73.3 



Table 7. Numerical values of minimum and maximum galacto- 
centric distances of the Sun under the action of galactic bar, disk 
+ halo and spiral arms. The four-armed spiral structure, v c = 
202.176 km s~ 1 , v y o = —220.0 km s~ 1 and v x o = 0.0 km s' 1 
are considered. The results hold for the time integration of 5 X 
10 9 years. The Sun remains (T / 100) X 5 X 10 9 years in the 
galactocentric distances Ft S (7.8, 8.2) kpc. 









Ftmax 


T 


[-] 


[ Vyo/Fto ] 


[ kpc ] 


[ kpc ] 


[%] 




U. i 


v 7ftn 
i . / oy 


Q Q A Q 


70 n 

( y.y 




U . o 


7 7ftO 

/ . / oy 


o.ooy 


en o 
oU.z 


n nn7 
U.UU / 


u.y 


7 77D 


Q QQQ 
O.OOO 


70 n 

/y.y 




1.0 


7.776 


8.340 


79.8 




1.1 


7.772 


8.333 


80.2 




0.7 


7.763 


8.372 


78.6 




0.8 


7.769 


8.366 


79.6 


0.014 


0.9 


7.766 


8.361 


79.1 




1.0 


7.778 


8.362 


78.3 




1.1 


7.767 


8.359 


78.9 




0.7 


7.758 


8.394 


77.7 




0.8 


7.763 


8.396 


77.9 


0.021 


0.9 


7.764 


8.389 


77.9 




1.0 


7.778 


8.385 


77.2 




1.1 


7.765 


8.376 


78.5 




0.7 


7.752 


8.422 


75.7 




0.8 


7.756 


8.419 


75.6 


0.028 


0.9 


7.759 


8.415 


76.1 




1.0 


7.771 


8.406 


76.2 




1.1 


7.762 


8.395 


77.2 



6 DISCUSSION 

The first case treated in the previous section, v c = 240 
km s _1 , corresponds to the value used by, e.g., Minchev 
& Famaey (2010), Famaey & Minchev (2010), Minchev et 
al. (2010), Minchev et al. (2011). The second approach char- 
acterized by v y o — —220.0 km s _1 is the conventional ap- 
proach, it uses the conventional speed (e.g., recommended 
by IAU) of the Sun 220.0 km s~ . As a consequence, the 
radial galactic migration of the Sun is smaller than ± 0.32 
kpc during the interval (5 — 10) x 10 9 years. Thus, the radial 
migration found by Minchev & Famaey (2010), Famaey & 
Minchev (2010), Minchev et al. (2010), Minchev et al. (2011) 
is caused not only by an incorrect equation of motion, but 
also by the fact that the effect of galactic bar together with 
other galactic components is sensitive to the real velocity of 
the Sun. 

Which of the two values of v y o is more realistic? We can 
answer this question by calculating the Oort constants. The 
first case, corresponding to v c — 240 km s , yields A — 17.1 
km s _1 kpc -1 and B — —14.9 km s _1 kpc' 1 , if the galactic 
disk, halo (Eq. I41[) and the dominant parts of the galactic 
bar (Eq. 01- $ (0) (r), $ (2) (r)) are considered. The second 
case, v y o = —220.0 km s _1 yields A = 15.0 km s" 1 kpc -1 



Table 8. Numerical values of minimum and maximum galacto- 
centric distances of the Sun under the action of galactic bar, disk 
+ halo and spiral arms. The four-armed spiral structure and v c 
= 202.176 km s _1 are considered. The results hold for the time 
integration of 5 X 10 9 years. The initial x— component of the 
Sun's velocity is v x . The negative sign of v x corresponds to the 
orientation toward the center of the Galaxy. The case when the 
resonance between the rotation of the Sun and spiral arms is con- 
sidered: e s = 0.028, Q s = 1-0 Vyo/Fto, v y o = 220.0 km s' 1 and 
Fto = 8 kpc. 





Rtmin 


Ftmax 


[km s' 1 ] 


[kpc] 


[kpc] 


+ 10.0 


7.558 


8.752 


0.0 


7.771 


8.405 


- 10.0 


7.774 


8.429 



and B = — 12.5 km s 1 kpc 1 . The first case is not consistent 
with the most probably values of the Oort constants, e.g., 
Oiling & Merrifield (1998) yield A = (11-15) km s" 1 kpc' 1 , 
B = -(12-14) km s' 1 kpc' 1 (see also Lequeux 2005, p. 10). 
Thus, the physical result should correspond to v y o = —220.0 
km s' 1 . Galactic migration of the Sun is less than ± 0.4 
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Table 9. Numerical values of minimum and maximum galacto- 
centric distances of the Sun under the action of galactic bar, disk 
+ halo and spiral arms. The two-armed spiral structure and v c 
= 202.176 km are considered. The results hold for the time 
integration of 5 X 10 9 years. The initial x— component of the 
Sun's velocity is v x . The negative sign of v x corresponds to the 
orientation toward the center of the Galaxy. The case when the 
resonance between the rotation of the Sun and spiral arms is con- 
sidered: e a = 0.06, f2 s = 1.0 Vyo/Ro, v y o = 220.0 km s _1 and Rq 
= 8 kpc. 





Rmin 


Rmax 


[km s^ 1 ] 


[kpc] 


[kpc] 


+ 10.0 


7.324 


8.436 


0.0 


7.333 


8.346 


- 10.0 


7.120 


8.667 



kpc with respect to the current value 8.0 kpc for the four- 
armed spiral structure. If we would add a nonzero value of 
radial velocity component to the velocity of the Sun, then 
the radial migration might be more relevant, see Table 8. 
However, the two-armed spiral structure yields much greater 
radial migration, up to 1.2 kpc for the largest value of the 
spiral structure strength. If we would add a nonzero value 
of radial velocity component to the velocity of the Sun, then 
the radial migration might be more relevant, see Table 9. 
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7 CONCLUSION 

The paper presents detailed derivation of the relevant equa- 
tions of motion for the Sun under the action of galactic bar 
and spiral arms. The equations differ from the equations 
used by other authors. 

Improving physics presented by Dehnen (2000), 
Minchev & Famaey (2010), Famaey & Minchev (2010), 
Minchev et al. (2010), Minchev et al. (2011) and others, 
we have obtained significantly different radial migration of 
the Sun than the previous authors. 

Galactic migration of the Sun is less than ± 0.4 kpc, if 
the four-armed spiral structure and the value v y o = —220.0 
km s _1 are considered. The Sun remains (69 — 79)% of 
its existence within the galactocentric distances (7.8 — 8.2) 
kpc and the result is practically insensitive on the spiral 
structure strength. Thus, the radial distance changes only 
within 5% from the value of 8 kpc. 

Galactic radial migration of the Sun is less than ± (0.3 
— 1.2) kpc for the two-armed spiral structure and v v o = 
-220.0 km s" 1 . The Sun remains (29 - 95)% of its exis- 
tence within galactocentric distances (7.8 — 8.2) kpc and 
the results strongly depend on the spiral structure strength 
and the angular speed of the spiral arms. The radial distance 
changes within (3.8 — 15.0)% from the value of 8 kpc. 



